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We have investigated the mechanism of stabilizing the simple-cubic (sc) structure in polonium 
(a-Po), based on the phonon dispersion calculations using the first-principles all-electron band 
method. We have demonstrated that the stable sc structure results from the suppression of the 
Peierls instability due to the strong spin-orbit coupling (SOC) in a-Po. We have also discussed the 
structural chirality realized in /3-Po, as a consequence of the phonon instability. Further, we have 
explored the possible superconductivity in a-Po, and predicted that it becomes a superconductor 
with T c ~ 4 K. The transverse soft phonon mode at q m |R, which is greatly influenced by the 
SOC, plays an important role both in the structural stability and the superconductivity in a-Po. 
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Polonium (Po), which belongs to chalcogen group in 
the periodic table, is unique in that it crystallizes in the 
simple-cubic (sc) structure [l| . Po having atomic number 
Z = 84 was discovered by Marie and Pierre Curie in 1898. 
Po exists in two metallic allotropes, a and /3-Po. It is a- 
Po that has a sc structure. Above 348 K, a-Po transforms 
to /3-Po, which has the trigonal structure [2|. Selenium 
(Se) and Tellurium (Te) , which are isoelectronic elements 
to Po in the chalcogen group, also have the trigonal spiral 
structure [3j. Note that the trigonal structure can be 
derived from the sc structure by elongation or contraction 
along the [111] direction [J], which is known to occur due 
to the Peierls distortion in p-bonded systems of Se and 
Te g. 

Recently, there have been several reports to explore the 
origin of the stabilized sc structure in Po jjMllj . General 
consensus so far is that the large relativistic effects in Po 
play an important role. The Peierls instability tends to 
be suppressed in Po by the relativistic effects. However, 
there was controversy on the role of spin-orbit coupling 
(SOC) of 6p electrons. Min et al.® claimed that the 
SOC in addition to the scalar-relativistic (SR) effects of 
mass-velocity and Darwin terms is essential, while Lcgut 
et aL@,0 claimed that the SR effects are already enough 
to stabilize sc-Po. The refined and comprehensive band 
calculations by the former group, by considering all the 
relativistic corrections and different exchange-correlation 
functionals show that the SOC of 6p states is really im- 
portant in stabilizing the sc phase of Po @ . 

In general, the structural transition is induced by the 
phonon softening instability, and so the phonon disper- 
sion calculations would provide evidence of the structural 
stability. Verstraeteflll calculated the phonon dispersion 
of sc-Po based on the pseudo-potential band method, and 
showed that the phonon softening instability does not oc- 
cur even without including the SOC. This result seems 
to support Legut et aVs claim[3] that the SR terms are 
sufficient to stabilize the sc phase. Thereafter, a couple 
of more studies were reported on the phonon dispersion 



of the a-Po at the ambient pressure [12j and under the 
pressure[l3|, by using the pseudo-potential band method. 
Since Po has a tiny energy difference, order of 1 meV, be- 
tween the sc and trigonal structures [8J, special caution is 
needed to calculate and analyze the phonon dispersions. 

The structural energetics in Po is known to be very sen- 
sitive to the volume and the utilized exchange-correlation 
method in the band calculation^, The lattice con- 
stant a exp of a-Po was first measured by Beamer et al. 14] 
to be 3.345 A. Later, Desando et a/.[H[ reported the 
refined lattice constant of a-Po to be 3.359 A. In the 
above band and phonon dispersion calculations, the old 
a °exp = 3.345 A was referenced. In the present study, we 
have taken a exp = 3.359 A as a reference. We have found 
that the small difference in a exp is important in analyzing 
the structural stability. 

As mentioned above, the Peierls instability is the key 
physics in Po ^EMl]. Even if the Peierls instability is 
suppressed by the SOC, a-Po in the sc phase would still 
have the soft phonon mode and so the strong electron- 
phonon (EP) coupling. Then the natural question is 
whether a-Po would have superconductivity or not. To 
our knowledge, superconductivity in Po has not been ex- 
plored yet, probably due to its radioactive and toxic na- 
ture. 

In this Letter, we have investigated phonon and super- 
conducting properties of a-Po, employing the first princi- 
ples all-electron band methods. The phonon dispersions 
and superconducting properties of a-Po are studied with 
and without the SOC at different volumes. By means 
of the phonon dispersion calculations, we have explic- 
itly demonstrated that the stability of the sc-phase arises 
from the SOC, which concludes the longstanding dispute 
on the origin of stable sc structure of Po. Further, we 
have predicted that a-Po is a superconductor with T c of 
- 4 K. 

We have employed the FLAPW band method 
implemented in the Elk package [13] ■ The SOC was 
included in sccond-variational scheme and the local- 
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FIG. 1: (Color Online) Phonon dispersions of sc-Po at a = 
a exp with SOC (LDA+SOC: blue-solid line) and without it 
(LDA: red-dotted line). The imaginary phonon frequencies 
are represented as negative values in the figure. 



density approximation (LDA) was used for the exchange- 
correlation energy. Actually, the choice of the LDA is 
more stringent test for the phonon stability than that 
of the generalized gradient approximation (GGA), be- 
cause the GGA always gives more softened phonons than 



the LDA. For sc-Po, a 



o _ 

th 



3.335A was obtained in the 
LDA+SOC scheme |l8j. which is in good agreement with 
o-exp — 3.359A of a-Po. Phonon dispersions were ob- 
tained by using the supercell method, implemented in 
the Elk package. The force constants and the dynam- 
ical matrix are obtained from the Hellmann-Feynman 
forces calculated with small individual displacements of 
nonequivalent atoms 

Figure [1] shows the phonon dispersions of a-Po at 
a — ci exp with and without the SOC [2(| ■ In the LDA 
scheme without the SOC (corresponding to the SR band 
scheme), the phonon softenings occur along all the high 
symmetry directions. Among those, phonon frequency 
along T-R (q » §R) is imaginary, which indicates the 



structural distortion along the [111] direction [2l|. In 
contrast, in the LDA+SOC scheme (corresponding to 
the fully- relativistic band scheme), the phonons become 
hardened and the imaginary phonon softening disap- 
pears. These features reflect that, without the 6p SOC, 
the sc structure of a-Po is unstable to a trigonal struc- 
ture, as in Se and Te. 

Note that, for the softened phonon at q w |R, there 
are two degenerate normal modes having opposite he- 
licities. The one has a clockwise helicity as shown in 
Fig. EJb), while the other has a counterclockwise helic- 
ity. Thus the softened phonon at q w |R induces the 
structural transformation to the trigonal /3-Po, which 
has the helical chain structure with definite chirality 
(Fig. HJc)). The chiral structure realized in /3-Po was 
once proposed in Te [22j, which was recently explained 
in terms of the orbital-ordered chiral charge-density wave 
(CDW) [H [13. This orbital-ordered CDW would be 
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FIG. 2: (Color Online) (a) One of two degenerate normal 
mode at q ~ |R in sc-Po. All the displacement vectors are 
placed on (111) planes, (b) The same normal mode shown 
along the [111] viewpoint, which manifests a clockwise he- 
licity. (c) The helical chain structure of /3-Po with definite 
chirality. Red and Green lines represent short bonds, (d) 
The difference of charge density of sc-Po in the LDA+SOC 
and that in the LDA (in unit of e/A ). Notice the charge 
density depletion in-between neighboring ions. 



Fig. [3] The chirality of similar feature was also observed 
in lT-TiSe 2 Interestingly, both Te and lT-TiSe 2 ex- 
hibit superconductivity under high pressure [26, 27 1 and 
Cu intercalation [28( , respectively, suggesting that the su- 
perconductivity emerges when the CDW is suppressed. 

The phonon softening in Fig. [T] is closely related to the 
Peierls mechanism. The phonon softening is described 
by the renormalization of phonon frequency by the EP 
interaction (the so-called Kohn anomaly), 



W 2 (q)=ft 2 (q)-|.g ep (q)| 2 X o(q), 



(1) 



suppressed, if the SOC were larger, as will be discussed in 



where w(q) and f2(q) correspond to renormalizcd and 
bare phonon frequencies, respectively, and g ep and Xo(q) 
are the EP coupling parameter and the charge suscep- 
tibility. Min et al. 6] obtained that Xo(q) m a-Po has 
the highest peak at q « |R due to the Fermi surface 
nesting. They also showed that the SOC suppresses the 
Fermi surface nesting effect and accordingly the Xo(q) 
value, whereby they predicted that the phonon softening 
is weakened and finally the imaginary phonon softening 
disappears. 

The suppression of the Peierls instability by the SOC 
occurs not only through the mechanism of Eq. (1) but 
also through weakening of the bonding strength. This 
feature is shown clearly in Fig. [2jd) , which plots the 
difference of the charge density in the LDA+SOC and 
that in the LDA. It is evident that, due to the SOC, 
the charge density is depleted in-between the neighbor- 
ing ions, which results in the weakening of the directional 
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FIG. 3: (Color Online) Phonon dispersion curves of sc-Te at 
the lattice constant of 3.210 A, without the SOC (LDA), with 
the SOC (LDA+SOC), and with the three times larger SOC 



bondings of Po chains. 

To verify the role of SOC in the structural stability 
more convincingly, we have examined the phonon dis- 
persions of Te in the hypothetical sc structure (sc-Te) 
with artificially varying the strength of SOC of Te (£so)- 
Figure [3] shows that, in the LDA, two imaginary phonon 
softenings occur along the T-M and T-R directions, while, 
in the LDA+SOC, that remains only along the T-R direc- 
tion. However, when we increases the SOC three times 
(3£|q), even the last imaginary phonon softening disap- 
pears eventually. This feature reveals that Te would also 
have a sc structure, if the SOC of Te becomes larger. The 
situation of the three times larger SOC in Fig. [3] actually 
simulates the case in Po, since the SOC of 6p electrons in 
Po (£so = 1-90 eV) is almost three times stronger than 
that of 5p electrons in Te (£,go = 0.72 eV). The behavior 
in Fig. |3] provide definite evidence that the SOC really 
plays a role of suppressing the Peierls instability. 

At the reduced volume, phonons tend to be hardened, 
and so it is expected that the phonon softening observed 
in Fig. [T] would be reduced. Indeed, the phonon disper- 
sions of sc-Po at a = 3.277 A (V/V exp = 0.93) in Fig.Hfa) 
do not exhibit the imaginary phonon softenings both in 
the LDA and LDA+SOC schemes. Figure 0|b) shows 
the behavior of phonon mode at q « f R as a func- 
tion of the lattice constant. It is seen that, with de- 
creasing the lattice constant, the phonon frequency rises 
very steeply in the vicinity of a exp (shaded area), and so 
the imaginary phonon softening disappears for a < a exp . 
This is the reason why the phonon softening instability 
does not occur at 



Hh 



3.335A (V/V exp = 0.98) of the 
LDA+SOC, as in the Verstraete's work [ll|, but occurs 
at a° th = 3.4111 (V/V exp = 1.05) of the GGA+SOC 0. 
Moreover, since the a to /3 transformation in Po occurs at 
348 K, one might have to take into account the thermal 
expansion (or) of the lattice constant [29l [ 3C | . Consider- 
ing ctT = 23.5 /im/m-K of a-Po at 298 K [3l|, the lattice 
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FIG. 4: (Color Online) (a) The phonon dispersion curves of 
the sc-Po with and without the SOC at the reduced volume 
V/Vexp = 0.93 (a = 3.277A). The filled circle corresponds 
to the phonon mode at q ~ |R. The overall phonon modes 
are hardened at the reduced volume, (b) The behavior of 
softened TA phonon mode at q « jR in the LDA with vary- 
ing the lattice constant. Here, a = 3.335A (V/V ex p = 0.98) 
corresponds to a^ h of the LDA+SOC, while a — 3.345A and 
a — 3.359A correspond to the low temperature experimen- 
tal lattice constants reported by Beamer et al. [1J| and by 
Desando et al. [l^|, respectively, a — 3. 366 A (filled square) 
corresponds to the lattice constant at 298 K [3ll |. 



elongation amounts to ~ 0.02 A at 300K, which is large 
enough to influence the phonon softening instability. In 
fact, the lattice constant of a-Po at 298 K is reported 
to be 3.366 A {V/V exp = 1.01) 0. Then, without the 
SOC, a-Po would be unstable to /3-Po, as indicated by 
the filled blue square in Fig. 0|b). 

We now consider the superconducting properties of the 
cv-Po [12] • In the Eliashberg theory, the average EP cou- 
pling constant is given by 



A = 2 / dtoa 2 F(u})/uj (2) 
Jo 

with the Eliashberg function expressed as 

a2F ^ = WE ) S l3kn,k+qm| 2 
^ ' qy k.n,m 

x <5(ekn)<5(ek+qm)£(w - Wqw). (3) 



Here N(E F ) is the DOS at E F 
frequency of the phonon mode v 
The EP matrix element g^ n k 



and Uqj, denotes the 
with momentum q. 
is obtained from 



qri 1 
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.9k„,k+qm = (kn|e q!/ -VV R rs(q)|k+q?n)/ v /2Ma; q „, where 
e qt/ is the polarization vector and VVifs(q) is the gra- 
dient of the Kohn-Sham potential with respect to the 
atomic displacements with the wave vector q. Then the 
critical temperature T c is obtained with the McMillan- 
Allen-Dynes formula [33| : 



LUl. 



°g 



1.20 



exp 



1.04(1 + A) 
A - n* - 0.62A^* 



(4) 



where wi og y = exp(| / °° ^a 2 F(w)lnw) j is the loga- 
rithmic average frequency and /i* is the effective Coulomb 
repulsion parameter. The estimated superconducting pa- 
rameters are listed in Table |TJ 

Noteworthy in Table Q] is that a-Po at the ambient 
pressure (V/V exp = 0.98) would be a superconductor 
with A = 0.89 and T c ~ 4K. At the reduced vol- 
ume (V/V exp = 0.93), both A and T c increase further 
to 1.04 and ~ 5 K, respectively. With reducing the vol- 
ume, N(Ep) decreases, while the Debye temperature 9d 
increases. The effect of the SOC on the superconduc- 
tivity is found to be not so large but detrimental. For 
V/V exp = 0.93, the SOC reduces A from 1.11 to 1.04, 
and T c from ~ 6 K to ~ 5 K. It is mainly because of 
the phonon hardening induced by the SOC, as discussed 
in Figs. Q] and HI This behavior in Po is quite opposite 
to that in neighboring element Pb, for which the SOC 
rather softens the phonons and increases A by more than 
40% (33 - l36j . This difference is thought to arise from the 
different Fermi surface nature between two 



34|. 



To examine which phonon contributes largely to su- 
perconductivity, we computed the q-dependent EP cou- 
pling constant A ql ,. As shown in Figs. [5ja)-(c), A qj ,'s 
show peaks at q's having phonon softening anomalies 
along T-X, M-R, and T-R directions. Longitudinal X 
and transverse T-R phonon modes have fairly large val- 
ues of A q „. We found above that A in sc-Po increases 
as the volume is reduced. In view of the phonon hard- 
ening at reduced volume, this behavior looks strange, at 
first glance. The enhanced A can be explained by the 
behavior of the Eliashberg function a 2 F(w) in Figs.[SJd) 
and (e). At the reduced volume, the maximum cut-off 
frequency is larger, and a 2 F(uj) is enhanced at low fre- 
quency. This leads to enhanced A and higher T c . As 



TABLE I: DOS at E F , N(E F ), Debye temperature 6 D , and 
the superconducting parameters of a-Po. Two values of T c 
are for two different values of /^*=0.10 and 0.13. 
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FIG. 5: (Color Online) The q-dependent EP coupling con- 
stants of a-Po: (a) in the LDA+SOC for V/V eq = 0.98, 
(b) in the LDA+SOC for V/V eq = 0.93, (c) in the LDA for 
V/Ve q = 0.93. The Eliashberg functions a 2 F(uj) with the 
SOC at different volumes of a-Po: (d) in the LDA+SOC for 
V/V eq = 0.98, (e) in the LDA+SOC for V/V eq = 0.93. Inte- 
grated Eliashberg functions A(w) are shown. 



shown in Figs. [SJd) and (e), the integrated Eliashberg 
function A(w) has a steep rise at about 4 and 5 meV for 
V/Vexp — 0.98 and V/V exp = 0.93, respectively. This en- 
ergy range corresponds to that of the transverse T-R soft 
phonon mode (Fig. 2]), suggesting that the superconduc- 
tivity of sc-Po comes largely from q « ^R phonon mode. 
In the LDA scheme of Fig. 03(c), the q s=s |R soft phonon 
mode produces the larger A qt/ , while, in the LDA+SOC 
scheme of Fig. EJb), the q « |R phonon mode hardened 
by the SOC reduces A q „ by half. This is the reason why 
the SOC is detrimental to the superconductivity in Po. 

In conclusion, based on the phonon dispersion calcula- 
tions, we have explicitly demonstrated that the SOC is 
the origin of the stabilized sc structure of a-Po. We have 
also predicted that a-Po would be a superconductor with 
T c ~ 4 K. The soft transverse phonon modes at q « |R 
in a-Po, which are greatly influenced by the SOC, play an 
important role both in the structural stability of the sc 
phase and its superconductivity. The experimental veri- 
fication of the superconductivity in a-Po is demanded. 
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